Phonon Spectra and Thermal Properties of Some fee Metals Using the 

Embedded- Atom Method 
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By employing the analytic embedded-atom potentials of Mei et al. [Phys. Rev. B 43, 4653 (1991)] 
we have calculated the phonon dispersion spectra for six fee metals: Cu, Ag, Au, Ni, Pd and Pt. We 
have also investigated thermal properties of these metals within the quasiharmonic approximation. 
Results for the lattice constants, coefficients of linear thermal expansion, isothermal and adiabatic 
bulk moduli, heat capacities at constant volume and constant pressure, Debye temperatures and 
Griineisen parameters as a function of temperature are presented. The computed results are com- 
pared with the available experimental data. The comparison shows a generally good agreement 
between the calculated and experimental values for all thermodynamic properties studied. Isother- 
mal and adiabatic bulk moduli and the specific heats are reproduced reasonably well, while the 
Griineisen parameter and Debye temperature are underestimated by about 10%. The calculated 
phonon frequencies for Ag and Cu agree well with the results from inelastic neutron scattering ex- 
periments. However, there is considerable room for improvement in the phonon frequencies for Ni, 
Pd, Pt and Au, particularly at high phonon wave vectors close to the Brillouin zone boundary. The 
coefficient of linear thermal expansion is underestimated in most cases except for Pt and Au. The 
results are good for Pt up to 1000K and for Au up to 500 K. 

PACS numbers: 63.20.-e, 65.40.-b 



I. INTRODUCTION 

The inadequacy of pure pair-potential models to de- 
scribe metallic cohesion is well-known and has been ad- 
equately documentedi^ Various approaches, at varying 
levels of sophistication, have been used to address the 
volume-dependence of the energy of a metallic system 
originating from the presence of the interacting electron 
gas. These range from using volume-dependent param- 
eters in the pair potential itself 3 to writing the total 
energy as a sum of pair potentials plus an empirical 
volume/density-dependent term££&Z& or an electronic 
band(bond) energy term^*A£ The latter is often written 
in terms of the first few moments of the electronic den- 
sity of states in the tight-binding approximationiii* 1 ^ For 
transition metals, invoking a simplified model of rectan- 
gular d-density of states due to Friedelji 3 - the bond energy 
term is sometimes approximated simply via the second 
moment of the density of statespi2*i£*i£*i£ 

An approach that has been widely used in this con- 
text is the embedded atom method (EAM) of Daw and 
Baskesji&iii where the energy of the metal is viewed as 
the energy to embed an atom into the local electron den- 
sity provided by the remaining atoms of the system, plus 
a sum of pair interaction potentials between the atoms. 
This method, which was developed almost simultane- 
ously by Finnis and Sinclair, 4 starts with the ansatz: 
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where E to t is the total internal energy, ph,i is the host 
electron density at atom i due to all other atoms, fj is 
the electron density of atom j as a function of distance 
from its center, i£y is the separation distance between 
atoms i and j. The function Fi(ph.i), called the em- 
bedding energy, is the energy to embed atom i in an 
electron density ph,i- &ij is a two-body central pair po- 
tential between atoms i and j. The host electron density 
Ph.i is assumed to be a linear superposition of spherically 
symmetric contributions from all individual atoms except 
atom i. 

In the model proposed by Finnis and Sinclair 4 the 
total energy of a system of atoms is assumed to con- 
sist of a binding term proportional to the square-root 
of the local density and a repulsive pairwise potential 
term. Another approach to this type of theory was pro- 
vided by Manninen, 20 Jacobsen, N^rskov, and Puska, 21 
who derived the functional form of the two terms in 
Eq.(l) by using the density functional theory. Although 
Daw and Baskes 1 ^*^ had also used the density func- 
tional theory to justify the use of Eq.Q, in most ap- 
plications of the method the two terms of Eq.Q were 
written in terms of parameters that were fitted to ob- 
served properties of the system. In this sense, the EAM 
remains an empirical or at best a semi-empirical method. 
However, the simplicity and ease with which it could 
be applied to a large variety of situations has led to 
its wide use in the study of liquidsj 22 i 23 i 24 alloyS) 25 i 26 
surfaces and interfaceS) 27 i 28 i 29 i 30 i 31 i 32 i 33 i 34 i 35 i 36 impurities 
and other defects in solids ULISLSL Relevant to the present 
work is the study by Foiles and Adams^ who investi- 
gated some thermodynamic properties of six fee metals: 
Cu, Ag, Au, Ni, Pd and Pt via quasiharmonic calcula- 
tions and molecular dynamics simulation using the EAM 
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potentials developed by Foiles, Baskes and Daw 
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II. THE MODEL 



By replacing the atomic electron density with an expo- 
nentially decaying function, JohnsorJi developed a set of 
analytic EAM functions for the nearest neighbor model of 
fee metals. However, this model had the limitation that 
all materials were forced to have the same anisotropy 
ratio of the shear moduli : 044/(011 — C12) = 1. Oh and 
Johnson^ extended the model beyond nearest neighours, 
at the cost of sacrificing the simple analytic form. Mei 
et al^ have overcome this challenge and extended this 
nearest-neighbor model into the one in which the embed- 
ded atom potentials are analytic and valid for any choice 
of the cut-off distance. In particular they derived the val- 
ues of the parameter for their EAM potential and density 
functions for six fee metals: Ni, Pd, Pt, Cu, Ag, and Au. 
They used this EAM potential in a molecular dynamics 
study of thermal expansion and specific heat of liquid Cu 
as a function of temperature, using a cut-off for the po- 
tential and embedding functions somewhere between the 
third and fourth neighbors for the corresponding crys- 
talline case. In later studies they applied the model to 
study self-diffusion in the liquid phase of the above six 
metals^ and the melting in AliSi Kuiying et a/4i have 
used this EAM potential in a molecular dynamics study 
of the local structure in supercooled liquid and solid Cu 
and Al. Although Mei et al~ had derived the values of 
the EAM parameters for six fee metals: Ni, Pd, Pt, Cu, 
Ag, and Au, to this date a systematic study of the vi- 
brational and thermodynamic properties of these solids 
in the bulk crystalline phase using these parameters has 
not been carried out. This work presents such a study 
for the first time and explores to what extent this partic- 
ular EAM model is successful in reproducing the exper- 
imentally observed phonon spectra and thermodynamic 
properties for the above six bulk fee metals. 



The organization of this paper is as follows: in section 
|n]we discuss some details of the EAM potentials of Mei et 
alm^ and its implementation in the present work. In sec- 
tion ^Q] we discuss the calculated phonon spectra for the 
six fee metals and their agreement with the experimen- 
tal results. In section llVl we discuss the calculated and 
the corresponding experimental values of various proper- 
ties: thermal expansion coefficients, isothermal and adia- 
batic bulk moduli, specific heats at constant volumes and 
pressures, the Debye temperature and the Griineiscn pa- 
rameter as a function of temperature. In section we 
summarize our results and conclusions about the validity 
of this particular EAM model. 



The embedding function F(p) and the two-body po- 
tential Q(R) in Eq. Q of Mei et al. are given by: 



F{p) = -E, 

+ ^2 s Aexp[-(p A - 1)7] 
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with 
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In Eq. ©, a is defined by 

a = (95 e O e /£ c ) 1/2 , 



(6) 
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£1 is the atomic volume, B the bulk modulus, E c the 
cohesive energy. The subscript e refers to the equilibrium 
value. s\ in Eq. JSJ is the number of atoms on the A-th 
neighbor shell with respect to a given reference atom. p\ 
refers to the A-th neighbor shell via 



Ra= PA Ri,A=1,2, 



(8) 



where R\ is the distance of the first neighbor shell with 
respect to a reference atom and R\ the distance of the 
Ath- neighbor shell. The constants pa depend on the crys- 
tal structure type: for the fee structure pa = VA- 

Mei et alm^ used n — 5 in Eq.© to fit the atomic 
charge density p. Only ratios of electronic densities ap- 
pear in Eq.© and hence p e cancels out. The constant f e 
was set equal to /0 e /12. The constant (3 was taken from 
previous works by Johnson^S^4 The remaining constants 
were obtained by fitting to the measured values of unre- 
laxed vacancy formation energy and the elastic constants, 
all of which were calculated by using three shells of neigh- 
bors for all the fee metals considered. The exact details 
of the fitting procedure used is given in Ref . ^(J . 

The values of the constants thus generated and given 
in Table I of Mei et aJ.^& were used in the present work 
with one notable difference. Mei et aJ4& obtained the val- 
ues of the parameters in their model by considering three 
shells of neighbors for the fee solid. In the molecular dy- 
namics study of liquid and solid Cu, they used a cut-off 



distance in their embedding function and pair-potential 
lying between the third and the fourth nearest neigh- 
bors. We have calculated the elastic constants Cu, C12 
and C44 using the homogeneous deformation method and 
also from the long wavelength phonons and studied their 
variation with respect to the number of neighbor shells 
included in the calculation. For three shells of neigh- 
bors our results obtained via the homogeneous deforma- 
tion method agrees well with those obtained by Mei et 
aJ.42i However, the values obtained via the homogeneous 
deformation method and from long wavelength phonons 
differ significantly from each other as well as from the 
experimental values. The differences in the results ob- 
tained by the two methods as well as between the calcu- 
lated and experimental values, decrease on increasing the 
number of neighbor shells in the calculation and practi- 
cally disappear as the number of neighbor shells reaches 
six. Based on this result, we have calculated the phonon 
spectra and all the physical properties by using six shells 
of neighbors for this model of EAM. We have compared 
the phonon frequencies for the test case of Cu by using 
three to six shells of neighbors The maximum difference, 
which is found for low energy phonons, between the three 
and six neighbor shells is of the order of 1%. The differ- 
ences in the phonon frequencies between the three and six 
neighbor shell calculations decrease with increasing wave 
vector, and at the zone boundary they are of the order 
of .01%. The static energy of the solid, given by Eq.JTJ, 
changes by ~ 0.7% as the number of shells increases from 
three to six. However, there is virtually no change in the 
location of the minimum in the static energy as a func- 
tion of lattice parameter, as well as the curvature at the 
location of the minimum. 
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FIG. 1: Phonon dispersion curves for Cu. The solid lines are 
the calculated phonon dispersion curves at the room temper- 
ature equilibrium lattice parameter of 3.6131 Ai— The square 
and round points are the experimental data at 296K from Ref. 
|4^ . L and T represent transverse and longitudinal modes, 
respectively. 



where I, m are the labels of the atoms. For E tot given by 
Eqs.(l) and (2), and with I ^ rn 
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III. PHONONS 



The phonon spectra were calculated as usual by diago- 
nalizing the dynamical matrix, obtained from the Fourier 
transform of the force constant tensor m) given by 



d 2 E t , 



dR l {l)dR3(m) 



(9) 



where the prime in the function in the above equation 
denotes the derivative of the function with respect to 
its argument. Though written somewhat differently the 
above expression is in agreement with that derived by 
Finnis and Sinclair^ 

The calculated phonon frequencies for the six fee met- 
als Ni, Pd, Pt and Cu, Ag, Au are plotted in Figs. lfTlrJjl . 
In Fig. Q we compare the calculated phonon spectrum 
in Cu with the experimental results of Svensson, Brock- 
house and flowed The squares and the circles represent 
the phonon frequencies from the inelastic neutron scat- 
tering experiment at room temperature (296 K), while 
the solid line denotes the calculated spectrum. For esti- 
mating the importance of the embedding term we have 
used dashed lines to show the frequencies obtained by 
considering only the pair potential term, i.e., neglect- 
ing the contribution from the embedding function. The 



4 



pair potential and the embedding terms contribute to the 
force constant with opposite signs: without the embed- 
ding term the phonon frequencies at high values of the 
wave vectors would have much worse agreement with the 
measured frequencies. We find that the three-body terms 
coming from the embedding function (the last term in 
Ea. (|10|l % l have a negligible effect on the phonon frequen- 
cies, the values obtained with and without these terms 
are virtually indistinguishable. We have also compared 
the calculated phonon frequencies with the experimental 
values for off-symmetry wave vectors, as given by Nilsson 
and Ronaldson.-i The agreement is similar to that for the 
wave vectors along the symmetry directions, in the sense 
that the agreement is very good for small wave vectors, 
while becoming progressively worse with increasing wave 
vectors to the same extent as for the symmetry direction 
wave vectors shown in Fig. Q). 

The term involving the embedding function in Eq.JQ) 
contributes about \ to | of the total energy for the model 
of Mei et al.^ but their contribution to vibrational fre- 
quencies is less than 10%. Hence, the embedding term 
contributes to the temperature variation of thermody- 
namic properties mainly via the static part of the free 
energy. 

For Cu the phonon frequencies calculated by consid- 
ering only the contributions from the nearest neighbors 
in the fee structure produce about 90% of the final con- 
verged frequencies originating from all neighbors. This 
is in agreement with the observation of Svensson et al.^ 
who used several Born-von Karman force constant mod- 
els involving neighbors up to various different shells to 
study the agreement between the calculated phonon fre- 
quencies and those from the inelastic neutron scattering 
experiment. Their least square fits of atomic and planar 
force constants to the observed phonon frequencies indi- 
cate that the nearest neighbor forces dominate in Cu, al- 
though longer range forces extending at least to the sixth- 
nearest neighbors are needed for a complete agreement 
between the calculated and experimental frequencies. A 
similar result was found by Nilsson and Ronaldsonii in 
a neutron crystal spectrometer study of phonon frequen- 
cies in Cu at wave vectors of both symmetry and off- 
symmetry points in the Brillouin zone at a temperature 
of 80 K. 

Among the six fee metals, the agreement between the 
calculated and experimental phonon spectra is almost 
perfect for Ag and only slightly less so for Cu. For all 
other metals the agreement becomes increasingly worse 
with increasing wave vectors in all the symmetry direc- 
tions, the worst case being that of Au. For Pd, Pt and 
Au the calculated phonon frequencies are underestimated 
with respect to the experimental results, while they are 
somewhat overestimated for Ni. For Cuj2& Ag^S Ni and 
Pd)2£ this trend as well as the level of agreement between 
the calculated and experimental phonon frequencies is 
the same as in the EAM scheme of Daw and Baskes^ 
where the total charge density at an atomic location is 
calculated from the ab initio Hartree-Fock results for free 
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FIG. 2: Phonon dispersion curves for Ag. The solid lines are 
the calculated phonon dispersion curves at the room temper- 
ature equilibrium lattice parameter^- The square and round 
points are the experimental data from Ref. |50(l at room tem- 
perature. L and T represent transverse modes and longitudi- 
nal modes respectively. 

atom charge densities4^& As pointed out by Daw and 
Hatcher ^22, as long as the fitting is done primarily to 
elastic constants, which involve phonons near the zone 
center only, the agreement for high frequency phonons 
near the zone boundary is not guaranteed. In this re- 
spect, phonons calculated via ab initio electronic struc- 
ture methods, which can capture the details of Fermi 
surface topology without resorting to any empirical fit- 
ting procedure, can yield superior results^2*i£ although 
results do vary depending on the method of electronic 
structure calculation and details of the implementation 
of exchange-correlation potentials, etc4^ Of course if one 
is simply interested in the phonon spectra, Born-von 
Karman fit to force constants can yield very satisfactory 
results4£ Note that the agreement with the experimental 
phonon spectra for Cu obtained by Cowley and Shukla 3 
by using a nearest neighbor Born-Mayer potential with 
volume-dependent prefactor is as good as obtained in the 
present EAM model. The total energy of the crystal in 
this study by Cowley and Shukla£ consisted of kinetic, 
exchange and correlation energies of the electron gas and 
an electron-ion interaction term in addition to the near- 
est neighbor Born-Mayer potential. 



IV. THERMODYNAMIC PROPERTIES 

In the quasiharmonic approximation, the total 
Helmholtz free energy of the crystal at temperature T 
and volume V or lattice constant a is given by 

F(a,T) = E tot (a) + F vib (a,T) 
= E tot (a) + k B 

Ks 

x]n{2sinh g-j, ' } (11) 
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FIG. 3: Phonon dispersion curves for Au. The solid lines, 
square and round points and the symbols L and T have the 
same meanings as in Figs. 1 and 2. The experimental data 
at 296K is from Ref. 10 



FIG. 6: Phonon dispersion curves for Pt. The solid lines, 
square and round points and the symbols L and T have the 
same meanings as in Figs. 1 and 2. The experimental data is 
from Ref. p| , 
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FIG. 4: Phonon dispersion curves for Ni. The solid lines, 
square and round points and the symbols L and T have the 
same meanings as in Figs. 1 and 2. The experimental data 
taken at 296K is from Ref. El. 
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where -Etot is the static total energy given by Eq. Q at a 
given volume V or lattice constant a. w s (K, a) is the fre- 
quency of sth mode for a given wave vector K and lattice 
constant a. U and Kb are the Planck's and Boltzmann's 
constants, respectively. 

All thermodynamic properties are calculated from the 
free energy given by Eq. . In obtaining the wave vec- 
tor sum in Ea. Hll(l we consider a 20 3 uniform grid in 
the Brillouin zone (BZ) giving 256 wave vectors in the 
irreducible part. The sum is performed using these 256 
wave vectors, appropriately weighted according to the 
point-group symmetry of the fee solid. Temperature vari- 
ation of thermodynamic properties is studied from OK to 
1400.fr: in steps of 2K. 

At a given temperature T, the equilibrium lattice pa- 
rameter is determined by the minmum of the Helmholtz 
free energy, i.e., 



,dF(a,T), _dE tot (a) ,dF vib (a,T), 



da 



da 



da 



= 0. (12) 



£ (K=[00Q2ji/a) £ (K=[0£l]2ji/a) £ (K=[0££]2ji/a) £ (K=[£££]2rc/a) 



The variation of lattice constant with temperature for 
the six fee metals is shown in Fig.J7J). It should be 
noted that our quasiharmonic results for Cu agree very 
well with those obtained by Mei et al~ using molecu- 
lar dynamics simulation. A small difference, increasing 
with temperaure, appears above 900 K due to the ana- 
harmonic effects not included in the quasiharmonic ap- 
proach. The results for all the six metals are as good as 
those obtained by Foiles and AdmasS using the EAM 
potentials of Foiles, Baskes and Daw 25 

The coefficient of linear thermal expansion is given by 



FIG. 5: Phonon dispersion curves for Pd. The solid lines, 
square and round points and the symbols L and T have the 
same meanings as in Figs. 1 and 2. The experimental data 
taken at 120K is from Ref. El. 
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In experimental works Eq. i|13fl is often replaced withS^ 
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FIG. 7: Lattice Constant Against Temperature for the fee 
Metals. Line a is for Ag, b is for Au, c is for Pt, d is for Pd, 
e is for Cu, and f is for Ni. Line g denotes the results of the 
MD simulation Mei et al. EH for Cu. 



a(T) 



da e (T), 



a e (T c ) dT 



(14) 



where T c is a reference temperature, usually taken to be 
the room temperature. To be consistent with experimen- 
tal and other theoretical work^ we take T c — 293. Using 
Eq. 114|) , we calculate the coefficients of linear expansion 
a(T) for the six fee metals and plot them in Figs.lpi- llU|) 
along with the experimental values. Thermal expansion 
coefficients for Cu, Ag and Ni calculated by Shukla and 
MacDonaldS^ and MacDonald and MacDonald^ based 
on empirical nearest neighbor central force model and 
incorporating two lowest order, cubic and quartic, an- 
harmonic effects show much better agreement with ex- 
periment than those given by the present EAM model in 
the quasiharmonic approximation. 

The isothermal bulk modulus at a temperature T is 
given by 



B(T) 



BIT) = V 



VI 



\ 9v )t 



dF 

p —w < 15 > 



dV 2 '"~ 
d 2 E to , 
dV 2 



V 



,d 2 F mb (a,T), 
dV 2 



(16) 



The second derivatives in Eq. IjlrJI) are obtained via nu- 
merical differentiation and the results are compared with 
the available experimental data in Figs. 1141 - I15fl . Both 
the values and the trend, a small decrease with increas- 
ing temperature, are reproduced reasonably well in the 
calculation. 
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FIG. 8: Coefficient of linear thermal expansion a(T) as a 
function of temperature for Cu. the solid line is the calcu- 
lated values, the square points are the experimental values 
from Ref. |57| . and the dashed line is the results of the MD 
simulation of Mei et al. liff 
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FIG. 9: Coefficient of linear thermal expansion a(T) as a 
function of temperature for Ag. The solid line denotes the 
calculated values, and th e sq uare points represent the exper- 
imental values from Ref. 
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FIG. 10: Coefficient of linear thermal expansion a(T) as a 
function of temperature for Au. The solid line denotes the 
calculated values, and th e sq uare points represent the exper- 
imental values from Ref. 15711 . 
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FIG. 11: Coefficient of linear thermal expansion a(T) as a 
function of temperature for Ni. The solid line denotes the 
calculated values, and the square points represent the exper- 
imental values from Ref. 15711 . 
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FIG. 12: Coefficient of linear thermal expansion a(T) as a 
function of temperature for Pd. The solid line denotes the 
calculated values, and th e sq uare points represent the exper- 
imental values from Ref. jo 71 . 
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FIG. 13: Coefficient of linear thermal expansion a(T) as a 
function of temperature for Pt. The solid line denotes the 
calculated values, and the square points represent the exper- 
imental values from Ref. 571. 
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FIG. 14: Isothermal bulk modulus B(T) as a function of tem- 
perature for Au, Ni and Pt. 
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FIG. 15: Isothermal bulk modulus B(T) as a function of tem- 
perature for Ag, Cu and Pd. 



The specific heat at constant volume Cy can be ob- 
tained from the temperature derivative of the total en- 
ergy. However, the contribution from the electrons ex- 
cited across the Fermi level cannot be obtained from the 
EAM expression. Hence we write the specific heat at 
constant volume as 



Cv{T) = C p v h {T) + C%{T). 



(17) 



where the phonon part is obtained from the EAM model 
via 



Cf(T) = ]TCV(K S ) 



K. 



Ks k 



2k B T 
1 



sinh 2 [hu s (K)/2k B T}' 



(18) 



We estimate the electronic part Cy(T) via the well- 
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known expression: 

CvP) = \k 2 B N{E F )T (19) 

where N{Ep) is the electronic density of states at the 
Fermi level Ep. In the above expressions we consider 
the specific heats and N(Ep) to be per atom. We 
use the well-known and well-documented Stuttgart TB- 
LMTO (Tight Binding Linear Muffin-tin Orbitals) codeS 
to compute the values of N(Ep) for various lattice pa- 
rameters. Our results agree well with those given by 
other methods. 59 With increasing temperature, lattice 
parameter increases and electronic bands become nar- 
rower, resulting in a slight increase in N(Ep). However, 
this increase is negligible and the variation of Cy (T) with 
temperature remains essentially linear. We have actually 
calculated the N(Ep) for 5-6 values of lattice parameters 
in the entire temperature range and used linear interpo- 
lation to find the values at all other lattice parameters. 
Ea. (|19fl . which is based on the independent fermionic 
quasiparticle picture^SiSi includes all the electron-ion, 
exchange and correlation effects as incorporated within 
the framework of density functional theory. What is ne- 
glected is the electron-phonon interactionjffiiSi which is 
presumably small for most of the metals considered, as 
none of these exhibits superconductivity down to almost 
absolute zero of temperature. Theoretical calculations 47 
of electron-phonon interaction show that for Cu it is in- 
deed negligible, while for Pd it may not be so. Since 
the values of the electron-phonon coupling constant are 
known only for a few of these metals, we have simply 
ignored this contribution. The experimental values con- 
tain the effects of electron-phonon interaction in addition 
to the contributions from vacancies 55 and other defects. 
Note that in the high temperature limit the calculated 
quasiharmonic value of C v h (T) reaches the classical har- 
monic value 3Ks per atom, as it should. 

We calculate the specific heat at constant pressure Cp 
by using the relation: 

C P (T)-Cy(T) = -T(^) 2 p (^) T . (20) 

or 

C P {T) = C v h (T)+C v l (T) 

+ ^a 2 (T)B(T)al(T)T. (21) 

Figs. (|16l - I21JI show the temperature-dependence of 
heat capacity at constant volume Cy , at constant pres- 
sure Cp and the electronic contribution Cy. The ex- 
perimental results, which are for Cp, are also shown. 
For Ni, Pd and Pt, where the Fermi level lies within (at 
the outer edge) the d-band, density of states N(Ep) is 
large, resulting in a substantial contribution to Cy- In 
most cases the agreement with the experimental results 
is good, except for Ni which is ferromagnetic below the 
Curie temperature T c of 627 Ki 4 ^ The discrepancy be- 
tween the calculated and experimental results is largest 
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FIG. 16: Calculated temperature-dependence of heat capacity 
of Ag at constant volume C v h , at constant pressure Cp and 
the electronic contribution Cy. The square points represent 
the experimental data for heat capacity at constant pressure 
from Ref. 0. 
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FIG. 17: Calculated temperature-dependence of heat capac- 
ity of Cu at constant volume C v , at constant pressure Cp 
and the electronic contribution Cy. The square points are 
the experimental data for heat capacity at constant pressure 
from Ref. |63j|. The dotted line denotes the results from the 
molecular dynamics simulation of Ref. |40| . 
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FIG. 18: Calculated temperature-dependence of heat capacity 
of Ni at constant volume C v h , at constant pressure Cp and 
the electronic contribution Cy. The square points are the 
experimental data for heat capacity at constant pressure from 
Ref. [it. 
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FIG. 19: Calculated temperature-dependence of heat capacity 
of Pd at constant volume C v , at constant pressure Cp and 
the electronic contribution Cy. The square points are the 
experimental data for heat capacity at constant pressure from 
Ref. El- 




FIG. 20: Calculated temperature-dependence of heat capacity 



of Pt at constant volume C v , at constant pressure C p and 
the electronic contribution Cy . The square points are the 
experimental data for heat capacity at constant pressure from 
Ref. El. 



U 2 



U 1 




FIG. 21: Calculated temperature-dependence of heat capacity 
of Au at constant volume Gy , at constant pressure Cp and 
the electronic contribution Cy. The square points are the 
experimental data for heat capacity at constant pressure from 
Ref. El- 
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FIG. 22: Adiabatic bulk moduli B(T) as a function of tem- 
perature for Ag, Cu and Pd. 



at the Curie temperature and decreases steadily for tem- 
peratures both below and above T c . For both low and 
high temperatures, away from T c , the agreement is very 
good. The discrepancy is understandable, as the present 
formulation of EAM does not distinguish between the 
magnetic and nonmagnetic states. For Cu (Fig. H17|) ~). 
the effects of anharmonicity can be seen in the dashed 
line, which represents the variation of Cp with tempera- 
ture obtained by Mei et al~ using molecular dynamics 
simulation. Note that despite poor agreement between 
the calculated and measured phonon spectra for Au, the 
calculated and measured values of Cp agree very well over 
a large temperature range, possibly due to anharmonic 
effects excluded from the present calculation, which are, 
however, manifested in all measured properties. Poor 
phonon frequencies and anharmonic effects compensate 
for each other, yielding excellent values for Cp. Although 
the calculated phonon spectra for both Au and Pt show 
poor agreement with measured values, measured Cp for 
Au shows much better agreement with the calculated val- 
ues than for Pt. On the contrary, measured thermal ex- 
pansion is much better reproduced for Pt than for Au by 
the calculation. 

The knowledge of Cp and Cy enables us to determine 
the adiabatic bulk modulus using the relation 



B S (T) = ^B(T). 



(22) 



Figs. il22l23|) show the calculated results. Experimental 
results are indicated wherever available. Both the values 
and the trend, a slight decrease with increasing temper- 
ature, are reproduced quite well by the model. 

An important thermodynamic property is the 
Griineisen parameter, which is essentially a measure of 
the volume-dependence of the phonon frequencies and 
can be related to thermal expansion^ of the solid. The 
mode-specific Griineisen Parameter corresponding to the 
(K, s) phonon mode is given by 



7s (K) 



V du a (K) 
"w.(K) dV 



(23) 
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FIG. 23: Adiabatic bulk moduli B(T) as a function of tem- 
perature for Au, Ni and Pt. 
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FIG. 24: Overall Griineisen Parameter 7 as a function of 
temperature for Ag, Cu and Pt. 




FIG. 25: Overall Griineisen Parameter 7 as a function of 
temperature for Au, Ni and Pd. 
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FIG. 26: Debye temperature On as a function of temperature 
T for Cu, Ag and Au. 



We compute the overall Griineisen parameter "f(T) by av- 
eraging over the individual Griineisen parameters 7 S (K) 
of all the modes with a weight of CV(Ks) from each 
mode(Ks)(2£ i.e., 



7(T) 



E Ks 7s(K)CV(K.s) 
£ Ks CV(Ks) 



(24) 



Figs. (|24l-l25p show the calculated results for the varia- 
tion of the overall Griineisen parameter j(T) as a func- 
tion of temperature. 

Finally, we study the Debye temperature and its 
temperature-dependence. The temperature dependent 
Debye temperature 0b is obtained by numerically eval- 
uating the integral 



C v (T)=9k B ( 



T > 



AT 4 exp(A:) 
[cxp(X) - !]■ 



■dX 



with 



x D = e D /T , 



(25) 



(26) 



and adjusting the value of Xjj so that the integral on 
the right hand side of Eq. I|25|) matches the already com- 
puted value of CV(T) on the left hand side. The results 



are plotted in Figs. (|2rj|> and i|27|). The calculated values 
are slightly smaller than the corresponding experimental 
values, reflecting the difference between the calculated 
and the measured values of Cy. However, the maximum 
difference is of the order of 10%. For all metals there is 
an initial decrease in the Debye temperature as the tem- 
perature increases from zero, except for Au where there 
is an initial increase. It is interesting that the present 
EAM model is able to capture this trend. 



V. SUMMARY OF RESULTS AND 
CONCLUSIONS 

This work presents the first complete study of the 
phonon spectra and all thermodynamic properties of six 
fee metals: Cu, Ag, Au, Ni, Pd and Pt, based on the 
EAM. The results are obtianed using the analytic em- 
bedding functions of Mei et ai4S in the quasiharmonic 
approximation. The calculated phonon dispersion curves 
for Cu and Ag agree well with the inelastic neutron 
diffraction results. The discrepancy between the cal- 
culated and the measured phonon frequencies increases 
with increasing phonon wave vector for all of the above 
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FIG. 27: Debye temperature <9_d as a function of temperature 
T for Ni, Pd and Pt. 

metals. However, the relative error, at the symmetry 
points X and L, is not more than 5.0% for Ag and 7.0% 
for Cu. Large differences between the calculated and the 
measured phonon frequencies are found for Ni, Pd, Pt 
and Au at high values of phonon frequencies, with the 
discrepancy for Au being the most drastic. However, the 
agreement or the lack thereof in the calculated phonon 
dispersion curves does not carry over to the thermody- 
namic properties in a proportionate way. For example, 
despite poor agreement for the phonon frequencies, the 
calculated and measured values of Cp agree very well for 
Au over a wide temperature range. Although the cal- 
culated phonon spectra for both Au and Pt show poor 
agreement with measured values, measured Cp for Au 
shows much better agreement with the calculated values 
than for Pt. On the contrary, measured thermal expan- 
sion is much better reproduced for Pt than for Au by the 
calculation. 

Both isothermal and adiabatic bulk moduli are very 
well represented over a wide temperature range for all the 
metals studied. This is perhaps not surprising, as the pa- 
rameters in the model are determined largely by fitting to 
the elastic constants, which usually have only very weak 
temperature-dependence. In general, the co-efficient of 



thermal expansion, calculated in the quasiharmonic ap- 
proximation, is underestimated for most of the metals 
beyond room temperature. Calculated values of Cp for 
Ag, Cu and Au agree very well with the measured val- 
ues. For Ni the agreement is good both below and above 
the Curie temperature, with the discrepancy increasing 
as the Curie temperature is approached from both above 
and below. This is understandable as the present model 
of EAM is not designed to capture the physics of a fer- 
romagnetic to paramagnetic transition. The Griineisen 
parameter and Debye temperature are underestimated by 
about 10% for all the metals. However, the temperature 
variation of the Debye temperature, an initial decrease 
as the temperature is raised from zero for all the met- 
als except Au, is reproduced well by the model. For Au 
the experimental values of the Debye temperature show 
an initial increase. This feature is also captured by the 
model. 

In summary, the thermodynamic properties of the six 
fee metals: Cu, Ag, Au, Ni, Pd and Pt, and their 
temperature-dependence are reproduced reasonably well 
by the EAM model of Mei et al.M The coefficient of ther- 
mal expansion is underestimated above the room tem- 
pearture in the quasiharmonic approximation used in this 
work. A better treatment of the anharmonicity in the 
potentials can improve the results, as evidenced by the 
molecular dynamics study of Mei et al.~ There is room 
for much improvement in the phonon dispersion curves 
for Ni, Pd, Pt and Au. There is some evidence that this 
can be achieved by a change in the fitting of the parame- 
ters. Preliminary calculations with small ad hoc changes 
in the parmeters representing the charge density in the 
model showed remarkable improvement in the phonon 
frequencies for Pd and Pt. This is an indication that a 
better representation of the charge desnity may lead to 
improvement for all of these metals. 
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